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-y ' Abstract 

This article discusses the properties of time evolution of wave pack- 
q-i ets in a few systems. Dynamics of wave packet motion for Rydberg 

atoms with the hierarchy of collapses and revivals is briefly reviewed. 
The main part of the paper focuses on the new mechanism of quan- 
tum recurrences in wave packet dynamics. This mechanism can oc- 
cur in any physical system with strong enough spin-orbit interaction. 
We discuss here the spin-orbit pendulum effect that consists in dif- 
ferent motions of subpackets with different spin fields and results in 
oscillations of a fraction of average angular momentum between spin 
and ordinary subspaces. The evolution of localized wave packets into 
toroidal objects and backwards (for other class of initial conditions) 
is also subject to discussion. 

PACS number(s): 03.65.Sq, 03.65.Ge, 32.90+a 
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1 Introduction 

Recent advances in experimental techniques have resulted in rapidly growing 
attention to the dynamics of particular nonstationary states, wave packets. 
With tunable, short and intensive laser pulses a creation and analysis of 
such states in atomic and molecular systems became available. An extensive 
research by means of both theoretical and experimental methods has brought 
an understanding of the intriguing phenomena of a hierarchy of collapses and 
revivals for particularly prepared wave packets in Rydberg atoms [1-7] and 
in Jaynes-Cummings model (JCM) of quantum optics [8-12]. 

In hydrogen and hydrogenoid atoms electrons can be excited to a coher- 
ent mixture of many Rydberg states that move almost classically for many 
Kepler periods [1-5,7]. However due to a nonequidistant spectrum of energy 
levels wave packets undergo a sequence of spreadings, collapses and revivals 
during long term evolution. There exist several time scales for the revivals 
(reccurences). The shortest one is associated with the period of a classical 
motion, i.e. Kepler period, the other ones called t rev and t sr (revival and su- 
perrevival times respectively) arise from subtle quantum interference effects. 
The examples of such effects are described in more detail in Section 2. 

In atomic systems the appropriate wave packet can be excited by a suit- 
ably chosen laser pulse [2,6]. The time evolution can be then probed by 
the following pulse after precisely measured time delay. Using additionally 
a circularily polarized electric field allows for a creation of electronics wave 
packets which are stationary [13,14]. In this case the wave packet motion is 
stabilized by a nonlinear coupling in the analogous way as orbits of Trojan 
asteroids are stabilized in the celestial mechanics [15]. 

Another interesting phenomenon already observed in quantum optics is 
an atomic analogue to a double-slit experiment [16]. Again in hydrogen or 
hydrogenoid atoms one can excite, by a pair of laser pulses with an appropri- 
ate time delay, two wave packets localized on the opposite sides of the orbit. 
Such a state is often called a 'Schrodinger cat state' as it is a mesoscopic size 
realization of the famous Schrodinger cat problem [17]. 

In the following sections 3 and 4 we discuss the new dynamical structures 
in wave packet motion. These new effects arise when the strong enough 
spin-orbit interaction is present in a physical system. In order to exhibit 
these effects in the most transparent way we have chosen the simple central 
potential in spherical harmonic oscillator form. This choice allows for the 



factorization of the evolution operator into two parts and for considering 
the orbital motion and spin-orbit motion separately, as well as for finding 
analytical solutions for both motions. In section 3 we present the spin-orbit 
pendulum effect [18-20] which is most pronounced for circular orbits. Section 
4 contains the discussion of linear trajectories for which spin-orbit interaction 
spreads (periodically) the initially spherical, well localized wave packet into 
a toroidal object, with a classical trajectory being the symmetry axis of the 
torus. Such behaviour of the wave packet has been presented here for the 
first time. 

All these experimental and theoretical achievements enrich our under- 
standing of the basic features of quantum mechanics. 

2 Properties of Rydberg wave packets 

The discussion presented in this section is in line of that of [3,5,7,22]. 

Consider a hydrogen atom prepared in such a way that the wave function 
is a superposition of bound-state eigenfunctions 

*(r, t) = J2 c n u n {r) exp {-iE n t/K) , (1) 

n 

where coefficients c n are not negligible only for rather narrow range of n 
around the mean value N. Obviously / u* n u s dr = 5 ns and J2n \ c n\ 2 = 1- This 
sum converges, therefore a finite number of c n makes it arbitrarily close to 
1 and is sufficient to represent \1/ with arbitrary accuracy. Then any ^ is 
arbitrarily close to a periodic function of time [23] . Let K be the common 
multiple of all n for which c n are not neglected. Then \1/ has a period T c iK 2 . 
The recurrence with this period is exact, but if many levels contribute to 
the sum (1) the required time is enormously long. However, nearly exact 
recurrences occur considerably earlier. The assumption that weights c n are 
strongly centered around a mean value N permits an expansion of the energy 
in a Taylor series 

E n ~E N + E' N (n-N)+ E% (n - N) 2 + E% (n - N) 3 + . . . , (2) 

where each prime on E^ denotes a derivative. The derivative terms in 
(2) define distinct time scales that depend on iV: T d = 2nN 3 , t rev = 
(N/3 + 1/2) T c i, t sr = (N 2 /4 + 1/2) T c i, called classical period, revival time 



and superrevival time, respectively. Keeping terms through third order and 
disregarding an overall time-dependent phase, one can rewrite (1) as 



*(r, t) =Y1 c nU n (r) exp 



J(n-N)t , {n-Nft t (n - N) 3 f 



■*■ cl ^rev v Sr J 

' (3) 

The convenient measure of the degree of recurrences is the recurrence 
probability (also referred to as autocorrelation function), where w n = \c n \ 2 
and time units are atomic units, 

p(t) = \(mmt))\ 2 = Y.™ne u / n2 . (4) 

n 

In order to exhibit some properties of short term and long term evolution let 
us construct in a standard way the so-called circular wave packet [1-4] . Then 
we can choose c n = (27T(t)~ 1//4 exp[— ((n— N)/2a) 2 } and u n (r) = $ n ,n-i,n-i(r) 
(aligned standard hydrogenic eigenfunction with I — m — n—1). The features 
of the evolution of such a circular wave packet are illustrated below with some 
numerical examples. We show the results obtained for the wave packet with 
N = 60, a = 1.5 and the summation over n taken from n\ = 50 to n,2 = 70. 
For this case t rev = 20.5T c i , t sr = 900.5 T c i and the mean principal quantum 
number N is large enough for the fractional revivals to appear. 

The short term evolution (time range 0-50 T c i) is displayed in Fig. 1, 
where the recurrence probability (4) is plotted as function of time. Initial 
peak intensity decreases very fast, and it is seen that after 4 revolutions 
wave packet is well spread over the whole orbit. Around t ~ 5T c i (1/4 of 
trev) peaks appear 4 times more frequently indicating formation of a frac- 
tional revival of 1/4 order (see also Fig. 3b). Around t ~ 7T C ; (1/3 of t rev ) 
and t ~ 10 T d (1/2 of t rev ) the frequency of peaks of the autocorrelation 
function indicate formations of the fractional revivals of orders 1/3 and 1/2, 
respectively (shapes of the wave packets are shown in Fig. 3c and 3d). The 
so-called full revivals are built at times close to t ~ 20 T c \ and t ~ 40 T c / and 
a corresponding shape of the wave packet is exhibited in Fig. 3e. 

Fig. 2 displays the same autocorrelation function (4) but in a much longer 
scale. The formation of recurrences corresponding to t sr at t ~ 300 i.e. 
t = t sr /3 and to some extent to t ~ 150 i.e. t = t sr /6 respectively is 
clearly seen. The corresponding shape of the wave packet at first superrevival 
occurence is present in Fig. 3f. 



3 Spin— orbit pendulum 

In this section we discuss the properties of the time evolution of wave packets 
representing a fermion moving in a spherical harmonic oscillator potential 
with a strong spin-orbit coupling. The hamiltonian of the system is then a 
simplified form of the single-particle Nilsson hamiltonian, extensively used 
in nuclear physics [24] 

H = H + V ls = H + K (T-a). (5) 

The extensive investigation of the problem is contained in refs. [18-21], here 
we review briefly the main results. In order to make the discussion simpler 
let us first focus attention on a particular case of circular orbits. Explicitly, 
we choose (without any loss of generality) the Oxy plane as the orbit plane 
and Ox as the initial spin direction (initial spin in orbit's plane). Then the 
initial states take the following explicit form 

Mt = 0)) = \N,a!) = \N)±(\+) + \-)), (6) 

where the eigenstate of s x is expressed explicitly by the eigenstates of s z 
operator (|+) and |— )) and \N) is the coherent state of spherical harmonic 
oscillator corresponding to a circular orbit. In configuration space it has the 
form 
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(f\N) = n~^ e~ 2l( x - x ° r+y + z J e ipw . (7) 

This packet has its maximum at (xo,0,0) and moves along the circular orbit 
if pq=xq=\N , N - the average value of n quantum number (in this case 
average values of l x , l y , l z are 0, 0, x^=N, respectively). Note that the length 
of the angular momentum is also N in these units. The state \N) can be 
decomposed in the HO basis as 

|JV)=£A,|ra = Z,Z,m = = £A,|m), (8) 

i i 

i.e. it is a linear combination of states with m—l and n—l. This particular 
state which combines all the partial waves and consequently evolves during 
the time evolution with the spin-orbit coupling by using all the frequencies of 
the spin-orbit partners, can be taken as 'pseudo' one-dimensional harmonic 



oscillator coherent state if one defines the weights A; in terms of I and a 
(continuous) real variable #asa Poisson distribution 



|Aj 



2 = e" N 



N l 

7T 



(9) 



As operators Hq and Vi. s commute, the evolution operator connected with 
the Hamiltonian (5) can be factorized as 



U(t)=U Q (t)U ls (t) 



e~ itH ° e ~ it( - l ' ff) 



(10) 



where the appropriate time units are chosen to absorb k and %. Using the 
equality (/ • a) 2 = I 2 — (I ■ a) in the expansion of Ui s (t) one can rearrange it 

as 

U ls (t)=e- it ^=f(t)+gm.a), (11) 

where, with the initial condition Ui s (t = 0) = 1, 

f(t) = s 

git) = e't 



cosiZ — sinXZ- 

2 n 2 



2i . t 



(12) 



(13) 



In the last equations Q=y/1 + 4/ 2 and states \lm) are its eigenstates with the 
eigenvalues (2/ + 1). 

The formulae (10-13) applied to (6) give the analytical form of the wave 
packet as a function of time. In the limit of high N (which in practice is 
achieved already for iV ~ 10) we can provide a simple explanation of the 
dynamics. In this limit we obtain \^(t)) = |\l/ + )|+) + |\l/_)|— ) where the 
components of the spinor are given as 



i*+> 



W 



-itH 



£A,e- ut |n/) 
i 



and 



1 ' V2 i 



(14) 



(15) 



Each component is again a coherent state that moves on a circular orbit of 
radius Xq, with angular velocities Uq — uji s and ujq + uji s , respectively. There is 
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in |^_) an overall phase e which produces at time t = -^Ti s = ir the simple 
result \^~) = — 1^+), which means that a pure state with spin along —Ox is 
then observed. 

These features are shown in details in Fig. 4 where the spatial motion 
of the wave packet (6-8) is plotted (more precisely only evolution under C/j s 
operator (11)). The figure displays the total probability density |\l/(t)| 2 for 
the case N = 18 at times £j = i/8 7] s , % — — 4. As motion is periodic the 
motion in the second half of the T\ s period is symmetric to that of the first 
half. Fig. 5 illustrates the motion in the spin subspace. It displays time de- 
pendence of the (s x ) and Tr(p 2 ), where p is the density matrix reduced to the 
spin subspace. The oscillations of (s) are accompanied by the correspond- 
ing oscillations of (I ) as the total angular momentum is conserved in the 
system. This effect of a periodic exchange of the portion of angular momen- 
tum between the spin and the ordinary subspaces has been called spin-orbit 
pendulum [18,19]. For analytical formulas and an extensive discussion see 
[18,19]. 

4 Time dependent vortex rings 

In this section we want to discuss the same physical system with the hamilto- 
nian (5) but quite different initial conditions corresponding to linear classical 
motion. More precisely we consider initial states in the form (6) but with 
the shape in configuration space at time t — given by 

(f\N) = 7T-2 e -5l(---o)W+^] _ ( 16 ) 

This is the gaussian wave packet at rest centered at (x ,0,0), x = v2N, 
which in HO potential moves along Ox axis with average value of the angular 
momentum equal zero (note that according to (6) spin direction is the same). 
The decomposition into the HO basis is more general 

\N} = J2^m m \nlm} . (17) 

nlm 

The coefficients X n i m are known analytically for arbitrary initial conditions 
(for details see eq. (55) of ref. [21] and related discussion). Again with (10- 
13) and (17) one can derive the analytical formulas for the time evolution. 



The details of this analytical formulation in terms of time dependent partial 
waves will be published elsewhere [25]. 

The examples of the motion of the wave packet with the lowest value 
of anglar momentum (L = 0) are shown in Figs. 6 and 7. Such a wave 
packet in absence of the spin-orbit coupling (k — 0) moves along the piece 
of the straight line on Oz axis. Both figures display only this part of the 
evolution which is governed by the Ui s operator, i.e. the harmonic oscillator 
motion is frozen. Both figures clearly show that the wave packet, initially 
well localized spreads into a torus (more precisely into two tori, the second 
with much lower density) with the classical trajectory being the symmetry 
axis. This spreading is reversible in our system due to its integrability and 
periodicity. Such 'vortex rings' topology of wave packets is presented here 
for the first time. The full motion is a superposition of spin-orbit motion 
and harmonic oscillator motion (10). 

5 Conclusions 

We have shown here that the presence of a spin-orbit interaction in a single- 
particle hamiltonian generates a new mechanism of collapses and revivals in 
wave packet dynamics. Assuming that initial state of a fermion is prepared 
in a form of a gaussian wave packet (a coherent state of harmonic oscillator) 
with well determined spin projection we are able to calculate its evolution 
in analytical terms. There exist some initial conditions for which formulas 
become compact and transparent, and can serve as a guidance for the most 
general cases. There are two extremes: the case N = L which corresponds 
to a circular classical orbit and the case L = corresponding to a linear 
orbit. In our simple model the constant form factor for spin-orbit coupling 
allows for a factorization of the evolution operator into two commuting parts. 
This fact makes possible considering both parts of the full motion separately 
and understand it better. The evolution operator Ui s , connected with the 
spin-orbit interaction results in splitting a wave packet into up and down 
components that evolve in a different manner. The so-called wave packet 
collapse reflects the vanishing overlap between subpackets. At (k + 1/2)TJ S 
the subpackets overlap again (exactly in the limit of high N), with the op- 
posite phases and the wave packet revives with spin reversed. At times kTi s 
the subpackets overlap exactly with the same phases and the initial state is 



exactly rebuilt. The sequence of collapses and revivals is accompanied by 
oscillations of average values of spin and orbital angular momentum, the so- 
called spin-orbit pendulum effect. This behaviour is common for any initial 
condition providing that spin is initially in orbit's plane. 

In the hydrogenoid atoms it is necessary to combine the properties of 
Rydberg wave packets described in section 2 with the spin-orbit evolution 
that we have derived. In this way we obtain a complicated interference patern 
in which the time scale due to the spin-orbit interaction must be included. 
No simple description has emerged yet. 
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Figure 3: 

Shapes of the wave packet (1) for t = 0,4.97,7.135,9.964,19.447,298.333 
(case a,b,c,d,e,f, respectively). Time instants chosen for this presentation 
correspond to sharp peaks of the autocorrelation function (4) displayed in 
Figs. 1 and 2. Therefore the part of the wave packet (for fractional revivals) 
or the whole packet (for full revivals) is localized in the same position as the 
initial wave packet. Coordinates are in units equal to the radius of Nth orbit. 
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Figure 4: Motion of the circular wave packet with N = 18 representing a 
fermion moving in the spherical harmonic oscillator potential with spin-orbit 
interaction (5). Shown is |*(t)| 2 = \^ + (t)\ 2 + \^-(t)\ 2 on the plane of the 
classical orbit (marked by the circle) for t{ = i T/ s , i=0-4 (from a) to ej) in the 
reference frame moving like the classical particle (in other words presented 
is the evolution due to Ui s only). 
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Figure 5: The time dependence of (s x ) and Tr(/9 2 ) (p - density matrix reduced 
to spin subspace). Time units equal to T/ s . 
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Figure 6: The same as in Fig. 4 but for initial conditions corresponding to the 
linear classical motion (iV=18, L=0) along Oz axis. Shown is the probability 
density of the wave packet on the Oyz plane. The full 3-dimensional wave 
packet possesses the cylindrical symmetry around Oz axis. HO evolution is 
frozen as in Fig 4. Cases a, b, c, d, e,f correspond to t 8 =0, 1/8, 2/8, 3/8, 7/16 and 
4/8T/ s , respectively. Note different vertical scales. 
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Figure 7: Cuts through the wave packet presented in Fig. 6 but by planes 
perpendicular to the symmetry axis with z = z cos t{ exhibiting a toroidal 
shape. Note different vertical scales. Time instants as in Fig. 6. 



16 



0.8 



0.6 



0.4 



0.2 




Figure 1: The recurrence probability (4) of the wave packet with N = 60, 
a = 1.5 as function of time (time unit equal T c \ - the period of a classical 
revolution). It is easy to recognize the structures of full revivals at t ~ 20, 40 
and fractional revivals at intermediate times. 
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Figure 2: The same as in Fig. 1, but for much longer times. Strong recur- 
rences at time t ~ 300 related to superrevival time t sr ~ 900 are clearly 
visible. 
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